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Abstract 

We study shear viscosities of different species in hot and neutrino-trapped dense matter relevant 
to protoneutron stars. It is found that the shear viscosities of neutrons, protons and electrons in 
neutrino-trapped matter are of the same orders of magnitude as the corresponding shear viscosities 
in neutrino-free matter. Above all, the shear viscosity due to neutrinos is higher by several orders 
of magnitude than that of other species in neutrino-trapped matter. 

Next we investigate the effect of shear viscosity in particular, neutrino shear viscosity on the 
thermal nucleation rate of droplets of antikaon condensed matter in protoneutron stars. The 
first-order phase transition from hadronic matter to antikaon condensed matter is driven by the 
thermal nucleation process. We compute the equation of state used for the calculation of shear 
viscosity and thermal nucleation time within the relativistic mean field model. Neutrino shear 
viscosity enhances the prefactor in the nucleation rate by several orders of magnitude compared 
with the T 4 approximation of earlier calculations. Consequently the thermal nucleation time in the 
T approximation overestimates our result. Furthermore, the thermal nucleation of an antikaon 
droplet might be possible in neutrino-trapped matter before neutrino diffusion takes place. 

PACS numbers: 97.60.Jd, 26.60.-c,52.25.Fi,64.60.Q- 
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I. INTRODUCTION 



A first order phase transition from nuclear matter to some exotic form of matter might be 
possible in (proto) neutron stars. It could be either a nuclear to quark matter transition or 
a first order pion/kaon condensation. Consequently, it might have tremendous implications 
for compact stars [l| and supernova explosions {2]. Here the focus is the first order phase 
transition proceeding through the thermal nucleation of a new phase in particular, antikaon 
condensed phase in hot and neutrino-trapped matter. After the pioneering work by Kaplan 
and Nelson on antikaon (K~ meson) condensation in dense baryonic matter formed in heavy 
ion collisions as well as in neutron stars J3I, several groups pursued the problem of antikaon 
condensation in (proto) neutron stars ^-18]. In most cases, the phase transition was studied 
using either Maxwell construction or Gibbs' rules for phase equilibrium coupled with global 
baryon number and charge conservation [19]]. The first order phase transition driven by 
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21] . In particular, 



nucleation of antikaon condensed phase was considered in a few cases 

the calculation of Ref. 2l| dealt with the role of shear viscosity on the the thermal nucleation 



of antikaon condensed phase in hot and neutrino-free compact stars [2l|. It is to be noted 



here that the first order phase transition through the therma 
droplets was also investigated in (proto)neutron stars 
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nucleation of quark matter 



nucleation theory of Langer 
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27] using the homogeneous 



2811 . The therma 



the quantum nucleation at high temperatures [25l. |27|. 



nucleation is an efficient process than 



We adopt the homogeneous nucleation theory of Langer [28|, [29] for the thermal nucle- 
ation of antikaon condensed phase. Nuclear matter would be metastable near the phase 
transition point due to sudden change in state variables. In this case thermal and quantum 
fluctuations are important. Droplets of antikaon condensed matter are formed because of 
thermal fluctuations in the metastable nuclear matter. Droplets of the new and stable phase 
which are bigger than a critical radius, will survive and grow. The transportation of latent 
heat from the surface of the droplet into the metastable phase favours a critical size droplet 
to grow further. This heat transportation could be achieved through the thermal dissipation 



and viscous damping 
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A parametrised form of the shear viscosity was used in earlier calculations of the nu- 
cleation of quark matter j^J. Recently, the influence of thermal conductivity and shear 
viscosity on the thermal nucleation time was studied in a first-order phase transition from 
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nuclear to antikaon condensed matter in hot neutron stars [21] . The shear viscosity due 
to neutrinos was not considered in that calculation. It would be worth studying the effect 
of shear viscosity on the thermal nucleation rate of droplets of antikaon condensed matter 
in neutrino-trapped matter relevant for protoneutron stars. Besides shear viscosities due 
to neutrons, protons and electrons, this involves the contribution of neutrinos to the total 
shear viscosity. Shear viscosities of pure neutron and neutron star matter were calculated 



by several groups 



32j-|37[]. We also performed the calculation of shear viscosity in neutron 



star matter using the equation of state (EoS) derived from relativistic field theoretical mod- 



els 



Transport properties of degenerate neutrinos in dense matter were estimated by 



Goodwin and Pethick 



We organise the paper in the following way. We describe models for EoS, shear viscosities 
of different species including neutrinos and the calculation of thermal nucleation rate in Sec. 
II. Results of this calculation are discussed in Sec. III. Sec. IV gives the summary and 
conclusions. 



II. FORMALISM 



The knowledge of the EoS for nuclear as well as antikaon condensed phases is essential 
for the computation of shear viscosity and thermal nucleation rate. We consider a first- 
order phase transition from the charge neutral and /3-equilibrated nuclear matter to K~ 
condensed matter in a protoneutron star. Those two phases are composed of neutrons, 
protons, electrons, electron type neutrinos and of K~ mesons only in the antikaon condensed 
phase. Both phases are governed by baryon number conservation and charge neutrality 



conditions 



19j . Relativistic field theoretical models are used to describe the EoS in nuclear 



and antikaon condensed phases. The baryon-baryon interaction is mediated by the exchange 



of a, u and p and given by the Lagrangian density 
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Cn = ipB (ilfid* 1 -m B + g aB cr - g^Bl^ - g P Bl^B • P M ) i>B 

B=n,p 

+ l - (d.ad^a - mla 2 ) - U(a) 
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The scalar self-interaction term ll|, |42|, |43) is 

1 



U(a) 



9i m N (g aN a) 



+ | 92 {g*Ncr) 



(2) 



The effective nucleon mass is given by m* B = itlb — gaB®, where mg is the vacuum baryon 
mass. The Lagrangian density for (anti)kaons in the minimal coupling scheme is given by 



c 



A 



ITKD^K - mfKK 



(3) 



where the covariant derivative is = + ig U K^^ + W p k^k ' P/j an d the effective mass of 
(anti)kaons is m* K = tjik — g<jK<3- For s-wave (p = 0) condensation, the in-medium energy 
of K~ mesons is given by 



u K - 



77 i 



K 



1 



(4) 



The condensation in neutrino-trapped matter sets in when the chemical potential of K~ 
mesons (hk- — ^k-) is equal to the difference between the electron chemical potential (/i e ) 
and neutrino chemical potential (fi Ve ) i.e. \Lk- — — Hu e - The critical droplet of antikaon 
condensed matter is in total phase equilibrium with the metastable nuclear matter. The 
mixed phase is governedby Gibbs' phase rules along with global baryon number conservation 
and charge neutrality 



We adopt the bulk approximation 22( which does not consider the variation of the meson 
fields with position inside the droplet. We solve equations of motion self-consistently in the 
mean field approximation (40| and find effective masses and Fermi momenta of baryons. 
Pressures in neutrino-trapped nuclear matter (P N ) and antikaon condensed matter (P K ) 
are given by Ref.ll3[. Here we calculate zero temperature equations of state because it was 
noted earlier that the temperature of a few tens of MeV did not modify the EoS considerably 
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Next we discuss the calculation of shear viscosity in protoneutron star matter. It was 
noted in earlier calculations that the main contributions to the total shear viscosity in 
neutron star matter came from electrons, the lightest charged particles, and neutrons, the 



most abundant particles 32 



34 l36l . l37f . Neutrinos are trapped in protoneutron stars and 



their contribution might be significant in transport coefficients such as shear viscosity. In 
principle, we may calculate shear viscosities for different particle species (n, p, e and v e ) 



in neutrino-trapped matter using coupled Boltzmann transport equations 
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38|. We 



can immediately write a set of relations between effective relaxation times (r) and collision 
frequencies (u^ ) following the Ref. 



2j 1 

i ,j=n,p,e,u e 



van + i^r,-) 



(5) 



which can be cast into a 4 x 4 matrix. However, solving the matrix equation becomes a 
problem because the relaxation time for neutrinos is much larger than those of other species 
as it is evident from Fig. 3. So we calculate shear viscosities due to neutrons, protons 
and electrons in neutrino-trapped matter from the above matrix equation as it was done in 
Ref. 38|. In this case, the effects of the exchange of transverse plasmons in the collisions of 
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44j. The knowledge o 



charged particles are included 
sections obtained from the Dirac-Brueckner approach [451 . 
of neutron and proton shear viscosities 



nucleon-nucleon scattering cross 
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46j is exploited for the calculation 



38l . 147] . On the other hand, we estimate the 



neutrino shear viscosity using the prescription of Goodwin and Pethick 
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k 2 (k + l) 2 [k(k + l)-2X v ] 



(6) 



For neutrino shear viscosity, we only consider scattering processes involving neutrinos and 



other species. Various quantities in Eq.(6) are explained in Ref. 39] and given below. The 
neutrino relaxation time (t u ) is, 



V r- 1 



i=n,p,e 



64tt 2 



< r > , 



(7) 
(8) 



and X v is defined as 



A, 



A! 



i 

J drt 2 dn 3 dn 4 \ [3(p\ ■ p 3 ) 2 - 1] (|M| 2 )^( Pl + P2 - P3 - P4 ) 
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10 /p 
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Here < |M 2 | > is the squared matrix element summed over final spins and averaged over 
initial spins for a scattering process and Cy and Ca are vector and axial vector coupling 
constants. 



For non-re 
reduces to 



39] 



ativistic nucleons {rrii/Ei) ~ 1, (pi/Ei) <C 1 and if also (p v /Ei) <C 1, \ l rj 



35 + 7°A t 



(15) 



However, we do not assume non-relativistic approximation in this calculation. The total 
shear viscosity is given by 

Vtotal = Vn + Vp + Ve + Vu , (16) 

where 

Vi(=n, P ,e) = -^r • (17) 

Effective mass 
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The relaxation time of i-th species (t^) is calculated using Ref. 
and Fermi momentum of i-th particle species are denoted by m* and pp v respectively. The 
electron effective mass is taken as its chemical potentials due to relativistic effects. 

We are interested in a first order phase transition driven by the nucleation of droplets 
of antikaon condensed phase in the neutrino-trapped nuclear matter. Droplets of antikaon 
condensed phase are born in the metastable nuclear matter due to thermal fluctuations. 
Droplets of antikaon condensed matter above a critical size (R c ) will grow and drive the 
phase transition. According to the homogeneous nucleation formalism of Langer and others, 
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the thermal nucleation per unit time per unit volume is given by 28|, |29 ] 

AF(R, 



r = r exp 



T 



(18) 



where AF is the free energy cost to produce a droplet with a critical size in the metastable 
nuclear matter. The free energy shift of the system as a result of the formation of a droplet 
is given by 

ITT „ v , (lg) 



AF(R) = - — (P K - P N )R 3 + AiraR 2 , 
3 

where R is the radius of the droplet, a is surface tension of the interface separating two 
phases and P N and P K are the pressure in neutrino-trapped nuclear and antikaon condensed 
phases, respectively as discussed above. We obtain the critical radius of the droplet from 
the maximum of AF(R) i.e. 5rAF = 0, 

2(7 



R 



c 



(P K - P N ) ' 

This relation also demonstrates the mechanical equilibrium between two phases. 



(20) 



We write t 



prefactors [29 



re prefactor in Eq. (18) as the product of two parts - statistical and dynamical 

3 

r = ■ (21) 

The available phase space around the saddle point at Rc during the passage of the droplet 
through it is given by the statistical prefactor (Oq), 
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2 © 3/2 (f)' 



(22) 



Here £ is the kaon correlation length which is considered to be the width of the interface 
between nuclear and antikaon condensed matter. The dynamical prefactor k is responsible 
for the initial exponential growth rate of a critical droplet and given by 
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2a 



K 



R 3 c (Aw) 2 



XT + 2(- V + C) 



(23) 



Here Aw = wk — is the enthalpy difference between two phases, A is the thermal 
conductivity and rj and ( are the shear and bulk viscosities of neutrino-trapped nuclear 
matter. We neglect the contribution of thermal conductivity because it is smaller compared 
with that of shear viscosity 2l|. We also do not consider the contribution of bulk viscosity 
in the prefactor in this calculation. 



We can now calculate the thermal nucleation time (r nuc ) in the interior of neutron stars 

as 

Tunc = {VVy 1 , (24) 

where the volume V = 4tt/3R^ uc . We assume that pressure and temperature are constant 
within this volume in the core. 



III. RESULTS AND DISCUSSION 



Nucleon-meson coupling constants in this calculation are taken from Glendenning and 
Moszkowski parameter set known as GM1 48j and those are obtained by reproducing the 
saturation properties of nuclear matter such as binding energy E/B = —16.3 MeV, baryon 
density n = 0.153 fm~ 3 , asymmetry energy coefficient a asy = 32.5 MeV, incompressibility 
K = 300 MeV and effective nucleon mass m* N /rriN = 0.70. 

Next kaon-meson vector coupling constants are estimated from the quark model and 



isospin counting rule 
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The scalar kaon-meson coupling constant is obtained from the 
real part of K~ optical potential depth at normal nuclear matter density. The strength 
of antikaon optical potential depth is obtained from heavy ion collision experiments and 



K atomic data 



49|-|54|. It has a wide range of values. On the one hand, the analysis 



of K atomic data predicted the real part of the antikaon optical pot ential to be as large 



as Uj{(no) = —180 ± 20 MeV at normal nuclear matter density Q, [50|. On the other 
hand, theoretical models including chirally motivated coupled channel models as well as a 
double pole structure of A(1405) yielded a less attractive antikaon optical potential depth 
Here we consider an antikaon optical potential depth of Ux(rio) = —120 MeV at 
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normal nuclear matter density and the corresponding kaon-scalar meson coupling constant 
is g a K — 1-6337. The value of antikaon optical potential adopted in this calculation resulted 
in a maximum neutron star mass of 2.08 M Q in earlier calculations using the Maxwell 



construction 



12] . This is consistent with the recently observed 2M & neutron star 591 ] . 



Using the above mentioned nucleon-meson and kaon-meson coupling constants, we calcu- 
late the EoS of neutrino-trapped nuclear and antikaon condensed phases at zero temperature 
in a self-consistent manner. A temperature of a few tens of MeV might modify the EoS at 
very high densities in (proto) neutron stars compared with the zero temperature EoS as it 
was noted in Ref . \v\ . However, the effect of finite temperature on the threshold of an- 
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tikaon condensation is negligible In this calculation, the EoS enters in Eq. (19) as 

the difference between pressures in two phases and in Eq. (23) as the enthalpy difference 
between two phases. Here we exploit the zero temperature EoS for the the calculation of 
shear viscosity and thermal nucleation time. The thermal nucleation of exotic phases was 



earlier investigated using zero temperature EoS in Ref. [2lL |26|. 

First we calculate shear viscosities of neutrons, protons and electrons in neutrino-trapped 
nuclear matter using Eq. (5) in the same fashion as it was done in Ref. 21J. We take lepton 
fraction Yl = 0.4 in this calculation. Shear viscosities of different species are shown as 
a function of normalised baryon density at a temperature T=10 MeV in Fig. 1. Here 
the electron viscosity is higher than the neutron and proton shear viscosities. The total 
shear viscosity excluding the viscosity due to neutrinos in neutrino-trapped nuclear matter 
is shown as a function of normalised baryon density for temperatures T = 1, 10, 30 and 100 
MeV in Fig. 2. The shear viscosity is found to increase with baryon density. Furthermore, 
the shear viscosity decreases with increasing temperature. The temperature dependence of 
shear viscosities is a complex one and quite different from the characteristic 1/T 2 behaviour 
of a Fermi liquid 36]. It is observed that the shear viscosities of neutrons, protons and 
electrons in neutrino-trapped nuclear matter are of the same orders of magnitude as those 
of the neutrino-free case 



Next we calculate the shear viscosity due to neutrinos. As a prelude to it, we compare ef- 
fective relaxation times corresponding to different species in neutrino-trapped nuclear matter 
in Fig. 3. Relaxation time is plotted with normalised baryon density at a temperature T=10 
MeV in Fig. 3. Relaxation times of different particle species due to scattering under strong 
and electromagnetic interactions are much much smaller than that of neutrinos undergoing 
scattering with other species through weak interactions. Consequently, particles excluding 
neutrinos come into thermal equilibrium quickly on the time scale of weak interactions. We 
calculate the shear viscosity due to neutrinos only treating others as background and it is 
shown as a function of normalised baryon density for temperatures T = 1, 10, 30 and 100 
in Fig. 4. Like Fig. 2, the neutrino shear viscosity decreases with increasing temperature. 
However, the neutrino shear viscosity is several orders of magnitude larger than shear vis- 
cosities of neutrons, protons and electrons shown in Fig. 2. It is the neutrino shear viscosity 
which dominates the total viscosity of Eq. (16) in neutrino-trapped matter. We perform 
the rest of our calculation using the neutrino viscosity in the following paragraphs. 
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We calculate the prefactor (To) according to Eqs. (21)-(23). The dynamical prefactor not 
only depends on the shear viscosity but also on the thermal conductivity and bulk viscosity. 
However, it was already noted that the thermal conductivity and bulk viscosity in neutrino- 
trapped nuclear matter were negligible compared with the shear viscosity 39|. We only 
consider the effect of shear viscosity on the prefactor. Besides transport coefficients, the 
prefactor in particular, the statistical prefactor is sensitive to the correlation length of kaons 
and surface tension. The correlation length is the thickness of the interface between nuclear 

re radius of a critical droplet is to 



and kaon phases [25|, |30j having a value ~ 5 fm [22J . T 



22 



301 ] . We perform our calculation 



be greater than the correlation length (£) for kaons 
with antikaon droplets with radii greater than 5 fm. The other important parameter in 
the prefactor is the surface tension. The surface tension between nuclear and kaon phases 



was already estimated by Christiansen and collaborators [60| and found to be sensitive to 
the EoS. We perform this calculation for a set of values of surface tension a = 15, 20, 25 
and 30 MeV fm -2 . The prefactor (To) is shown as a function of temperature in Fig. 5. It 
is shown for a baryon density n\, = 4.235no which is just above the critical density 3.9rio 
for antikaon condensation at zero temperature [l3], and surface tension o = 15 MeV fm~ 2 . 
The prefactor was also approximated by T 4 according to the dimensional analysis in many 
calculations [26|, [30j. The upper curve in Fig. 5 shows the prefactor of Eq. (21) including 
only the contribution of neutrino shear viscosity whereas the prefactor approximated by T 4 
corresponds to the lower curve. It is evident from Figure 5 that the approximated prefactor 
is very small compared with our result. 

Now we discuss the nucleation time of a critical droplet of antikaon condensed phase in 
neutrino-trapped nuclear matter and the effect of neutrino shear viscosity on it. The thermal 
nucleation rate of the critical droplet is calculated within a volume with R nuc = 100 m in 
the core of a neutron star where the density, pressure and temperature are constant. The 
thermal nucleation time is plotted with temperature for a baryon density = 4.235no in 
Fig. 6. Furthermore, this calculation is done with the kaon correlation length £ = 5 fm and 
surface tension o = 15, 20, 25 and 30 MeV fm -2 . The size of the critical droplet increases 
with increasing surface tension. Radii of the critical droplets are 7.1, 9.4, 11.7 and 14.1 
fm corresponding to a = 15, 20, 25 and 30 MeV fm -2 , respectively, at a baryon density 
4.235no- The nucleation time of the critical droplet diminishes as temperature increases for 
all cases studied here. However, the temperature corresponding to a particular nucleation 
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time for example 10~ 3 s, increases as the surface tension increases. There is a possibility that 
the condensate might melt down if the temperature is higher than the critical temperature. 
So far, there is no calculation of critical temperature of antikaon condensation in neutrino- 
trapped matter. However, the critical temperature of antikaon condensation was investigated 
in neutrino- free matter in Ref. 18) . We compare thermal nucleation times corresponding to 
different values of the surface tension with the early post bounce time scale td ~ 100 ms in 
the core collapse supernova [2] when the central density might reach the threshold density 
of antikaon condensation. The time scale td is much less than the neutrino diffusion time 
~ 1 s as obtained by Ref. [3^]. Thermal nucleation of the antikaon condensed phase may be 
possible when the thermal nucleation time is less than td- For a = 15 MeV fm -2 , the thermal 
nucleation time of 10~ 3 s occurs at a temperature 16 MeV. It is evident from Fig. 6 that 
the thermal nucleation time is strongly dependent on the surface tension. Further thermal 
nucleation of an antikaon droplet is possible so long as the condensate might survive the melt 
down at high temperatures [18]. Our results of thermal nucleation time are compared with 
the calculation taking into account the prefactor approximated by T 4 in Fig. 7 for surface 
tension a = 15 MeV fm -2 and at a density rib = 4.235no- The upper curve denotes the 
calculation with T 4 approximation whereas the lower curve corresponds to the influence of 
neutrino shear viscosity on the thermal nucleation time. The results of the T 4 approximation 
overestimate our results hugely. For a nucleation time of 10~ 3 s at a temperature T=16 Mev, 
the corresponding time in the T 4 approximation is larger by several orders of magnitude. 



IV. SUMMARY AND CONCLUSIONS 



We have studied shear viscosities of different particle species in neutrino-trapped (3- 
equilibrated and charge neutral nuclear matter. We have derived equations of state of 
nuclear and antikaon condensed phases in the relativistic mean field model for the calcula- 
tion of shear viscosity. It is noted that neutrons, protons and electrons come into thermal 
equilibrium in the weak interaction time scale. The shear viscosity due to neutrinos is 
calculated treating other particles as background and found to dominate the total shear 
viscosity. 

Next we have investigated the first-order phase transition from neutrino-trapped nuclear 
matter to antikaon condensed matter through the thermal nucleation of a critical droplet of 
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antikaon condensed matter using the same relativistic EoS as discussed above. Our emphasis 
in this calculation is the role of the shear viscosity due to neutrinos in the prefactor and its 
consequences on the thermal nucleation rate. We have observed that the thermal nucleation 
of a critical antikaon droplet might be possible well before the neutrino diffusion takes place. 
Furthermore, a comparison of our results with that of the calculation of thermal nucleation 
time in the T 4 approximation shows that the latter overestimates our results of thermal 
nucleation time computed with the prefactor including the neutrino shear viscosity. Though 
we have performed this calculation with antikaon optical potential depth of U^(no) = —120 
MeV, we expect qualitatively same results for other values of the antikaon optical potential 
depth [61 1 . 
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Fig. 1. Shear viscosities corresponding to different particle species in neutrino-trapped nuclear 
matter are shown as a function of normalised baryon density at a temperature T = 10 MeV and 
Y L = 0.4. 
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FIG. 2. Total shear viscosity in neutrino-trapped nuclear matter except the contribution of 
neutrinos is plotted with normalised baryon density for different temperatures. 
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Fig. 3. Relaxation times corresponding to different species in neutrino-trapped nuclear matter are 
shown as a function of normalised baryon density at a temperature T = 10 MeV and Yj, = 0.4. 
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Fig. 4. Neutrino shear viscosity is shown as a function of normalised baryon density at different 
temperatures. 
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FIG. 5. Prefactor including the contribution of shear viscosity is plotted as a function of temper- 
ature at a fixed baryon density and surface tension and compared with that of T 4 approximation. 
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FIG. 6. Thermal nucleation time is displayed with temperature for different values of surface 
tension. 
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FIG. 7. Same as Fig.6 but our results for a fixed surface tension are compared with the calculation 
of T 4 approximation. 
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